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Abstract. We observed SN 1979C in M100 on 4 June 1999, about twenty years after explosion, with a very sensitive 
four-antenna VLBI array at the wavelength of A18cm. The distance to M100 and the expansion velocities are 
such that the supernova cannot be fully resolved by our Earth-wide array. Model-dependent sizes for the source 
have been determined and compared with previous results. We conclude that the supernova shock was initially in 
free expansion for 6±2 yrs and then experienced a very strong deceleration. The onset of deceleration took place 
a few years before the abrupt trend change in the integrated radio flux density curves. We estimate the shocked 
swept-up mass to be M sw ~ 1.6Mq, assuming a standard density profile for the CSM. Such a swept-up mass for 
SN 1979C suggests a mass of the hydrogen- rich envelope ejected at explosion no larger than M onv ~ 0.9Mq. If 
SN 1979C originated in a binary star the low value of M env suggests that the companion of the progenitor star 
stripped off most of the hydrogen-rich envelope mass of the presupernova star prior to the explosion. 
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1. Introduction 

Supernova SN1979C in M100 was disc overe d on 1979 
April 19 by Gus E. Johnson (Mattei et al. |l979| ) at around 
magnitude 12 in the visible, and it is t hough t to have 
exploded on April 4, 1979 (Weiler et al. |l986| ). Optical 
spectra first showed a featureless continuum, whic h late r 
evolved to exhibit strong Ha emission (e.g., Schlegcl |l996j ). 
SN1979C was extraor dinari ly luminous, M^ ax w -20 
(e.g., Young & Branch |1989| ), making it among the most 
luminous type II supernova ever observed. SN 1979C was 
also found to be a radio supernova (RSN) (Weiler and 
Sramek, 1980). Indeed, at a distan ce of 16.1 ± 1.3 Mpc 



estimated for M100 (Ferrarese et al. 1996) it is one of the 
intrinsically strongest RSN (Weiler et al. 1996). From Ha 



widths an expansion speed of 9 200±500kms was esti- 



mated (Panagia et al. 1980 and G. Vettolani, priv. comm.) 
for epoch around 45 d. 
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SN1979C has been classified as a SN of type II-L. 
Fesen et al. (1999), hereafter F99, have studied the late- 
time optical properties of SN 1979C and other three type 
II-L supernovae. Late-time optical emission lines from 
type SNe II-L appear surprisingly steady over long time in- 
tervals. These authors find similar expansion velocities in 
the range 5100-6000 km s _1 in the strongest lines present 
(Ha, [Oi] 6300, 6364 A, [On] 7319, 7330 A, and [Om] 
4959, 5007 A) for all type II-L supernovae. From the asym- 
metries of the lines these authors find evidence of dust 
formation. Strong radio emission from a SN appears to be 
strongly correlated with late time bright optical emission, 
with SN 1979C as a prime example. In SN 1979C, all de- 
tected lines but Ha show emission peaks at —5 000 and 
— 1000 km s -1 suggestive of clumpy emission from mate- 
rial coming mostly from ejecta, and physically separated 
from the material emitting Ha. The Ha emission tracks 
well the peak 6 cm emission (Weiler et al. 1991 ) indicating 
that both emissions are related to the material from the 
swept-up shell. 
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Fig. 1. Radio light curves for SN 1979C in M100 at wave- 
lengths of 20, 6, and 2 cm. The curves represent the best- 
fit model light curves. The best-fit parameters were de- 
termined using only data through 1990 December (day 
~4300). Notice the flux density increase from day 4000 
onwards (Montes et al. 2000). 



The evolution of the radio emission from SN 1979C has 
been interpreted within the mini-shell model (Chevalier 



1982 ) as due to synchrotron emission from the outer part 



of the shock swept-up shell of circumstellar material at- 
tenuated by thermal absorptio n from the even mo re dis - 
tant ionized CSM (Weiler et al. |1986| , Montes et al. |2000D . 
Synchrotron self absorption which has been show n to be 
relevant in SN1993J (Fransson & Bjornsson 1995, Perez- 
Torres et al. |2001|) migh t also play a role in SN 1979C, 
although Chevalier ( 1998 ) is not of this opinion. The mod- 
ulation present i n the radio light curves o f SN 1 979C has 
led Weiler et al. Jl992| ) and Montes et al. ( fcOOOD to model 
the progenitor as possibly being in a detached eccentric 
binary system with a less massive companion. ^From a 
study of the emission environment of SN 1979C, Van Dyk 
et al. (1999) have estimated the mass of the progenitor 
to be 17 — 18 (±3) M©. Hydro-dynamical simulations by 
Schwarz & Pringle ( 1996| ) confirm that a binary system 
similar to that proposed by Weiler et al. (1992) with a 
spiral shaped stellar wind around the progenitor is feasi- 
ble for SN 1979C. Such modulation of the density of the 
gas surrounding the RSN and now being shocked, might 
have dramatic effects on the shape of the radio structure. 
A determination of such radio structure would be a pow- 
erful way to further understand the pre-supernova phase 
of the progenitor of SN 1979C. 

The radio emission of this RSN, monitored by Montes 
et al. ( pOOC ), has started to show an interesting trend 
which can be seen in Fig. |l| after years of steady de- 
cline the radio emission has been constant, or perhaps 
increasing, since about 1990 (~ 10 yr) in a way somewhat 
reminiscent of the radio emission increase for the case of 



SN 1987A (Ball et al. |1995| ). It appears that after a stan- 
dard RSN radio emission decline phase (Weiler et al. 1996) 
the expanding shock may have encountered denser ionized 
gas which has boosted the radio emission. These results 
may be in accord with those of F9£ which show that high 



electron densities (10 cm -3 ) are necessary to explain 
the absence of several lines, such as [On] 3726, 3729 A in 
their optical spectra. 

Soft X-rays have been detected from SN 1979C by 



Immler et al. ( [1998D , Ray et al. ( p001[ ), and Kaaret fl200ip . 
The soft X-ray luminosity appears constant or slightly de- 
clining from epoch 16 yr through 20 yr. This result appears 
compatible with the model predic tions of circumstellar in- 
teractions (Chevalier & Fransson 1994 ). 

VLBI observations at 3.6, 6, 13 and 18 cm by Bartel et 
al. (1985), hereafter B85, determined a size of SN 1979C 
and its growth for epochs ranging from Dec 1982 to 
May 1984 assuming a synchrotron emitting optically thin 
spherical model. At those early epochs the source could 
not be resolved by any VLBI array. At an average ex- 
pansion speed of, for example, ^8 000 kms -1 , one should 
expect a shell radius growth rate of ~ 0.3/iasday" 1 which 
over 20 years should yield a shell radius of ~2.2 mas which 
could be resolved by A6 cm VLBI if detected on the long 
baselines. However, due to sensitivity limitations of the 
available arrays, the chances of detecting it at A 6 cm at 
the present time are null. On the other hand, at A18cm 
the emission strength of SN 1979C and the sensitivity of 
the arrays make the detection possible. We have carried 
out state-of-the-art A18cm observations in an attempt to 
learn as much as possible about the source structure. 

2. Observations 

SN1979C was observed on 4 June 1999 from 19:30 UT 
to 01:30 UT at 1.6 GHz with a sensitive transatlantic 
VLBI array. The antennas used were (diameter, loca- 
tion): Effelsberg (100m, Germany), Robledo (70m, Spain), 
phased- VLA (130m-equivalent, NM, US), and Goldstone 
(70m, CA, US). Data were recorded using a standard 
VLBA mode with a data rate of 256Mb/s, covering 
a bandwidth of 64 MHz. The data were correlated at 
the VLBA Correlator Center of the National Radio 
Astronomy Observatory (NRAO)[}n Socorro, NM, US. 

The observations were made with a duty cycle of about 
26 min consisting of 3 min observations of the calibrator 
source PKS B1157+156, followed by 5 min observations of 
the calibrator source TXS 1214+161, in turn followed by 
14 min observations of the source SN1979C. The source 
3C274 was observed twice as a calibrator, once at the 
beginning and once at the end of the observations. 

The correlated data were analyzed using the 
Astronomical Image Processing System (aips). A ba- 

1 NRAO is operated under cooperative agreement with the 
National Science Foundation. 
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sic integration time of 2 sec was used. The phase slopes 
across each of the 8 MHz bands and over the 64 MHz 
synthesized band were aligned using a phase calibration 
determined from observations of the calibrator source 
3C274. The data of the calibrators PKSB1157+156 
and TXS 1214+161 were then fringe-fitted in a standard 
manner. In the next step, phases, delays, and delay- 
rates corresponding to the sources PKS Bl 157+156 and 
TXS 1214+161 were interpolated to the source SN 1979C, 
to provide an initial model for the fringing of the data 
of this radio source. A solution was obtained by fring- 
ing the data using narrow search windows set around 
the initial model. Fringe detections of SN 1979C were ob- 
tained for all baselines at almost all times. However, some 
SN 1979C data had to be removed from the analysis after 
visual inspection of the data of the two main calibrators. 
Indeed, for some time ranges in the observations there 
were calibration problems or data with lower than accept- 
able quality. For instance, high winds at the phased- VLA 
produced an auto-stowing of some antennas and problems 
in the phasing during the second half of the observations. 
Therefore, data from the VLA from 23:30 UT to 01:15 UT 
were edited. 

Adjacent to the VLBI observations reported in this 
paper, some cross-scans on PKSB1157+156 were carried 
out at 1.6 GHz to fix the pointing parameters of the an- 
tenna in Effelsberg. As a byproduct, those observations 
allowed the total flux density of PKSB1157+156 to be 
measured as 0.79+0.02 Jy (A. Kraus, priv.comm.). This 
value was adopted to set our flux density scale. Then, by 
readjusting the gains of the antennas, the VLBI map of 
PKSB1157+156 was forced to have the same flux density 
as the total flux density (that is, PKSB1157+156 was as- 
sumed unresolved by our VLBI array.) With this calibra- 
tion procedure for the antennas of the array a total flux 
density of 5.4+0.2 mJy was obtained for SN1979C. 

Further model fitting analysis was carried out with 



modelfit in difmap (Shepherd et al. 1994), and with 



QFIT in the Caltech VLBI Analysis Package (Pearson & 



Readhead 1984). 



then many models which can reproduce the data equally 
well (see Marscher |1985| and |B85| ). 

Thus, we decided to determine the size of SN 1979C 
for the following putative source shapes: an optically thick 
source, an optically thin shell of width 30% of the exter- 
nal radius, and a ring that is an infinitely bright shell 
of zero width. The last one is obviously unphysical and 
its determination is only valuable as a limit. Even the 
first of the three putative source shapes is not expected 
given that the spectrum of the source is optically thin with 
S cx ;,-0-63±o.03 (Montes et al. |2000D- The diameters deter- 



mined by model fitting from our data for such models are 
4.57 + 0.25, 3.60 + 0.17, and 3.10 ± 0.14mas, respectively. 
The uncertainties were approximately estimated from the 
statistical errors given by the programs QFIT and mod- 
elfit, taking into account discrepant results yielded by 
them. 

A fit for the preferred model of an optically thin shell of 
diameter 3.60 ± 0.17mas with width 30% of the external 
radius is shown, along with the measured visibilities, as 
the solid line in Fig. |^. 



SNI979C at I ■'.Hi 999 Jun 04 



Fig. 2. Visibility amplitude versus interferometric dis- 
tance. The model (solid line) of 3.60 mas outer diameter 
shell with a 30% thickness is plotted superimposed on the 
data. 



3. The size of SN 1979C 



We assumed spherical symmetry in the modeling. The 
data remaining after the editing were of the highest qual- 
ity, but from a look at the data it was not fully clear that 
all points in the visibility vs. uv-radius plot are on the 
main lobe of the interferometric response. We had to prove 
to ourselves that such was indeed the case. Good measure- 
ments of the closure phases (all being around zero and 
none around 180°) were a great help. We conducted sev- 
eral tests and model fits to convince ourselves that indeed 
the array at A18 cm was not large enough to be able to re- 
solve the structure of SN 1979C. Unfortunately, there are 



4. Discussion 

Marcaide et al. 
SN1993J a shell 



(1995a ( 


1995b 




1997) 



determined for 
structure with a shell width 30% of 
the outer radius, somewhat larger than that predicted by 



Chevalier & Fransson (1994). Therefore, we have adopted 
a shell of width 30% of the external radius as a model 
for SN 1979C in this paper and all our discussion will be 
based on such model. For the present discussion all pre- 



vious results by B85 will be transformed (normalized) to 
this model. 
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At a distance of 16.1 Mpc to M100 an angular size of 
1 mas corresponds to a linear size of 2.41 • 10 17 cm. Hence, 
the source outer radius (1.80 mas) corresponds to 4.33 • 
10 17 cm. Since 20.12 yr have elapsed since the explosion, 
the corresponding average supernova expansion velocity 
is ^6 800 km s^ 1 . Since early expansion speeds were esti- 
mated at epoch 45 d from the Ha lines as ^9 200 km s _1 it 
follows that the expansion has been strongly decelerated 
by the action of the circumstellar matter. How strongly? 
It is difficult to make a detailed analysis to answer such a 
question, since intermediate size estimates are not avail- 
able. Nonetheless, we now discuss our VLBI results, in 
conjunction with previously published VLBI and optical 
spectroscopy data, to show that SN 1979C has indeed suf- 
fered a heavy deceleration in its expansion, starting 6 ± 2 
yr after its explosion. 

B85 determined for SN 1979C a diameter of 
2-16iQ'5Qinas 5.17yr after explosion, from 1.67 GHz 
VLBI observations using an optically thin uniform 
sphere model. This value translates into a diameter of 



1.56 



+0.31 
-0.36 



mas for our discussion shell model. This size 



corresponds to an average speed for the first 5.17yr 
of ^11500kms _1 . On the other hand, the average 
speed of the supernova over the first six weeks after 
explosion as determined from the widths of Ha lines was 
- 9 200 ± 500 kms" 1 (Panagia et al. |1980| , G. Vettolani, 
priv. comm.), which is about 20% lower than the average 
velocity estimated for the first 5 years from the B85 
estimates normalized to our discussion model. These two 
results are nicely consistent. Indeed, in the mini-shell 
model of Chevalier 1982, the radio emission comes from 
the forward shock, while the optical, UV, and X-ray 
emissions come mainly from the reverse shock, which 
moves slower than the forward shock. Hence the velocity 
ratio determined from the VLBI and optical observations 
is roughly what is expected for the radii ratio of the 
forward and reverse shocks, in a model of 20-30% shell 



width (Chevalier & Fransson |1994| , Marcaide et al. [1997D 
if SN 1979C followed a free, self-similar expansion for 
the first 5 years after its explosion. A free expansion for 
the first 4 years, followed by a decelerated expansion 
would also be compatible with the mentioned results. 
These results should be understood as evidence of free 
expansion of SN 1979C for the first 4-5 years and of the 
presence of a shell-like structure, and not as evidence of a 
given shell width, since different model normalizations of 
the results of B85 would yield equally consistent results 
and a similar conclusion. 

An estimate of m (R oc t m ) using the size estimate of 



B85 and our size estimate results in m 



0.62+°;^. This 
results shows that SN 1979C has indeed dramatically de- 
creased its growth rate after a likely initial phase of free ex- 
pansion (to sa 1.0; see Fig. ||). The optical data, although 
as scarce as the VLBI data, can be used as a check of the 
reliability of our estimate of to. In a self-similar expan- 
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Fig. 3. Angular diameter vs. age of SN 1979C. B85 in- 
dicates the size estimate from Bartel et al. (1985) nor- 
malized to our model (see text), F99, the size estimated 
by us based on the optical results by Fesen et al. (1999), 
and M02 indicates our results. The solid line indicates 
a possible expansion, which is free for the first 5 yrs and 
decelerated from then on. The size estimate F99 shown 
is compatible with this expansion scenario and it would 
be slightly different for scenarios where the free expansion 
would last much longer than 5 yr, which is not the case. 
See text. 



sion scenario R oc t m and, consequently, v oc i m_1 . If we 
now assume that the supernova did not decelerate for the 
first 5.17 years, we can take the value of 9 200 kms -1 as 
the expansion velocity just prior to the phase of strongly 
decelerated expansion. Combining this value with that of 
6 200±300 kms" 1 at epoch 14.09 yr (F99 and R. Fesen, 
priv. comm.), gives to = in agreement with our 

estimate of m = 0.62^'^- 

Hence, m = 0.62 appears to be an adequate value to 
characterize a phase of strong deceleration in the expan- 
sion of SN1979C, assuming it started about 5.17 years 
after its explosion. If the supernova started to decelerate 
later than 5.17yr, then to should be smaller, and corre- 
spondingly stronger the deceleration of SN 1979C. It seems 
somewhat arbitrary to assume that the deceleration set in 
at the same time as the observations of B85. Indeed, there 
is nothing apparent in the integrated radio light curves 
which could point to a change at that particular epoch. It 
is intriguing, though, that the integrated radio light curves 
of SN 1979C (Fig. 1) seem to show an increase in flux at 
around t=2200 days, or 6 yr after its explosion. One could 
speculate that the flux rise at this epoch is somehow re- 
lated to a change in the deceleration parameter, likely due 
to the shock front of the supernova entering a high-density 
region of the circumstellar medium (CSM). 

Our estimate of the angular diameter for SN 1979C 
can also be used to obtain an upper limit for the epoch 
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at which the phase of strong deceleration started. We find 
that the latest epoch of free expansion which is still com- 
patible with both the optical observations of F99 (epoch 
14.09yr) and our VLBI determination at 20.12 yr, and 
with the fact that a lower limit is set by SN 1979C enter- 
ing a Sedov phase (m = 0.4), corresponds to ~8 yr. 

The integrated radio emission curves have been well 
modeled for the first 10 years by Weiler et al. ( 1992| ) and 
Montes et al. (200C) who assumed self-similar expansion 



and a modulation in the density of the CSM due to the 
binary interaction, to explain the m odulation in the flux 
density curve. Montes et al. ( |2000|) found that starting 
around epoch 10 yr the radio emission started to grow af- 
ter a long decline predicted by the radio emission model. 
Can our VLBI observations shed any light o n wh at hap- 
pened at epoch 10 yr? The VLBI results by B85 and us 
require a strong deceleration of the expansion. Such de- 
celeration is compatible with the optical results by F99 
for a self-similar expansion. Our results rule out any ini- 
tial free expansion longer than 8 years. Thus, the onset 
of decelerated expansion took place well before the epoch 
10 yr where the radio emission decline turned onto a slight 
growth. The sooner the onset of the decelerated expansion 
the milder the deceleration needed to be compatible with 
the observations. The epoch of that onset is restricted to 
the range 4-8 years and the corresponding deceleration 
parameter m to the range 0.67-0.44. (Strictly speaking an 
onset at 3 years with a corresponding m=0.72 is not ex- 
cluded, but it appears unlikely since it yields a radio to 
optical velocity ratio of 0.68 at epoch 5.17 yr. Such a ratio 
is small, but not yet inconsistent with the models used. 
In any case, the deceleration required for an early onset 
is comparable to the stronges t dece lerations reported in 
other cases by Marcaid e et al . (|l997|) for SN 1993J in M81 
and Mc Donald et al. ( |200l| ) for 43.31+592 in M82). An 
onset of decelerated expansion at epoch 6 yr (correspond- 
ing to a deceleration parameter m=0.57) appears most 
likely and it may be associated with a small halt in the 
otherwise decreasing flux density trend, a halt which may 
be partly related to encountering a denser CSM. 

If the expansion of SN 1979C has significantly deceler- 
ated, as our VLBI estimate shows, then the mass of the 
CSM swept up by the shock front, Af sw , must be compa- 
rable to or larger than the mass of ejected hydrogen-rich 
envelope, Af onv . Our estimated angular diameter trans- 
lates into a forward shock radius of ~ 4.33 x 10 17 cm. 
For a presupernova wind of velocity i> w =10kms -1 , this 
radius corresponds to ^13 700yr prior to the explosion of 
SN 1979C. Weiler et al. (1991) estimated a mass-loss rate 
of - 1.2 x 10~ 4 M Q vr" 1 for the progenitor of SN 1979C. 
If we now assume that the circumstellar medium up to 



this distance is a standard one (p c 



oc r 



~ 2 ), we then ob- 



tain AL, 



1.6Mq. Momentum conservation therefore 



implies that M c 



0.59Af sw , which results in a mass for 



the ejected hydrogen-rich envelope of ~ 0.94Af ( 



©• 



In the case of the Type II supernova SN1993J, whose 
progenitor was a ~ 15Mq star in a binary system, mass 
exchange with its companion meant that the mass of 
the hydrogen-rich envelope ejected at the explosion was 
0.2 - 0.4M© (Woosley et al. 1994, Houck & Fransson 
1996), much less than the w 3.3A/q one would expect 
from a single-star model (Hoflich et al. 1993). Weiler et 
al. (1986) have suggested that SN1979C was born in a 
binary system consisting of a 15-18 Mq red supergiant 
and a 10 Mq Bl main-sequence star. If we assume that 
SN 1979C lost, like SN 1993J, most of its hydrogen-rich 
envelope prior to explosion, then the ~ 0.9Mq deceler- 
ated ejecta suggested from our observations represents at 
least all of the remaining hydrogen-rich envelope. 

/From VLBI observations of SN 1979C more than 20 
years after its explosion, we conclude that an initial free 
expansion of SN 1979C was followed by a strong decelera- 
tion which set on ~6 yr after the supernova explosion. We 
estimate the shocked swept- up mass to be M sw ~ 1.6Mq, 
assuming a standard density profile for the CSM. Such 
a swept-up mass for SN 1979C suggests a mass of the 
hydrogen-rich envelope ejected at explosion as low as 
A/cnv ~ 0.94Mq or lower. If, a s sugg ested by Weiler et al. 
flL986D and Schwarz & Pringle ( |1996| ) SN 1979C originated 
in a binary star (alike SN 1993J) this value of Af env sug- 
gests that the companion of the progenitor star stripped 
off most of the hydrogen-rich envelope mass of the pre- 
supernova star prior to its explosion. New observations 
at A18 cm around 2005 should be able to test the above 
conclusion and would also be essential to image the true 
structure of SN 1979C. On the other hand, a continued rise 
in the level of flux density at A6 cm and improved VLBI in- 
strumentation and techniques like phase-referencing may 
allow an earlier determination of the intrinsic structure of 
SN 1979C with high resolution and a test of the binary 
star scenario which the present results tend to favor. 
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